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SECTION  I 


INTRODUCTION 


The  explosive  material  used  in  the  largest  quantity  per  year  in  the  world  is  ammonium  nitrate 
(AN)  mixed  with  a  near  stoichiometric  amount  (about  6  wt%)  of  hydrocarbon  fuel  oil  (FO). 
Millions  of  pounds  of  this  mixture,  denoted  ANFO,  are  produced  each  year  and  used  for 
about  80  percent  of  all  commercial  blasting  operations.  The  US  consumption  of  ammonium 
nitrate  based  explosives  was  4.0  billion  pounds  in  1991,  accounting  for  97  percent  of  all 
industrial  explosives  that  were  used.  Large  quantities  of  ANFO  are  therefore  handled  in 
manufacturing  and  transportation  operations.  The  thermal  stability  and  initiation  insensitivity 
of  ANFO  are  very  high,  making  ANFO  much  safer  to  handle  and  store,  especially  when 
compared  with  previously  used  explosives  such  as  dynamite.  However,  this  high  degree  of 
safety  may  have  contributed  to  a  relaxation  of  handling  precautions — possibly  too  large  a 
relaxation. 

The  purpose  of  this  study  was  to  quantify  the  thermal  hazards  inherent  in  using  large  charges 
of  ANFO  in  various  DNA  testing  programs.  For  example,  over  the  past  several  years  ANFO 
charges  from  about  2400  to  4800  tons  (2.2  -  4.4  x  1 06  kg)  have  been  used  to  simulate  the 
environments  caused  by  nuclear  weapons.  In  hemispherical  shapes,  these  weights  represent 
quantities  of  ANFO  with  radii  of  1070  and  1350  cm  (70  and  88  foot  diameters).  The  main 
question  addressed  in  this  work  was  whether  it  is  possible  for  ANFO  charges  of  this  size  to 
self-heat  as  a  result  of  chemical  decomposition  heat  which  cannot  be  conducted  away  as  fast 
as  it  is  generated. 
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SECTION  2 

BACKGROUND 

2. 1  THERMAL  DECOMPOSITION  KINETICS  FROM  COOK-OFF  DATA. 

Thermal  methods  of  measuring  reaction  rates  rely  on  the  heat  produced  (or  absorbed)  by  a  system  as  a 
measure  of  the  amount  of  reaction.  We  have  utilized  thermal  runaway  (isothermal  cook-off)  tests  on 
ANFO  to  determine  critical  temperatures  (i.e.,  minimum  temperatures  for  runaway  self-heating)  from 
which  we  have  derived  steady-state  global  reaction  rates.  Since  different  sample  sizes  yield  different 
critical  temperatures  (i.e.,  different  rates  of  reaction),  we  have  been  able  to  determine  the  global  rale 
constant  parameters  for  ANFO  decomposition  over  a  range  of  over  live  decades. 

Thermal  runaway  from  exothermic  reactions  is  possible  for  many  chemicals  and  energetic  materials. 
Tliis  happens  when  the  rate  of  heat  generation  from  chemical  decomposition  exceeds  the  rate  of  heat 
loss  from  the  material.  At  the  critical  temperature,  these  two  rates  are  equal.  Quantitatively 
describing  the  rate  of  heat  loss  by  conduedon  from  simple  sample  configurations  is  a  relatively 
straightforward  problem;  once  the  problem  is  solved,  the  rate  of  heat  generadon  can  be  determined. 

Since  the  rate  of  heat  generadon  is  the  product  of  the  reacdon  rate  and  the  heat  of  reacdon,  the 
reacdon  rate  can  be  determined  if  the  heat  of  reacdon  is  known.  Cook-off  tests  on  samples  of 
approximately  25  kg  or  so  yield  rate  constant  deterrainadons  in  the  lO'V  range,  whereas  tests  on 
smaller  samples  yield  rate  data  at  much  higher  values,  about  HE6  and  10 :  s'1  for  1  kg  and  40  mg 
samples,  respectively.  The  corresponding  rate  constant  at  the  cridcal  temperature  for  approximately 
2,000  tons  would  be  in  the  10  “  s'1  range. 


2.2  THERMAL  EXPLOSION  MODELS. 


The  theory  of  thermal  ignition  has  been  treated  by  many  workers,  including  Zeldovitch,  et  al.1  (1939), 
Frank- Kamenetskir  (1939),  Gray  and  Lee3  (1968),  Zinn  and  Mader4  (1960),  Zinn  and  Rogers5  (1962), 
and  many  others.  For  solid  explosives,  consideradon  can  be  confined  to  conducdve  heat  losses  [i.e., 
the  well-known  Frank- Kamenetskii  (F-K)  modeil.  The  assumptions  and  limitadons  of  this  model  have 
been  discussed  by  the  authors  listed  above.  In  the  case  of  symmetric  conilguradons,  the  F-K  formula 
gives  the  cridcal  temperature  (Tc,  K)  as 

Tc  =  Ea/R  , 
ln\r±pQAEa\ 

Tc2  X  ShR 


where 


Critical  temperature,  K 
Activadon  energy,  cal/mol 
Prcexponenlial  factor,  s'1 
Gas  constant,  1.987  cal/K/mol 
Radius  or  half-thickness,  cm 
Density,  g/errr 


2 


Q  =  Heat  of  decomposition  reaction,  cal/g 
k  =  Thermal  conductivity,  cal/cm/s/K 

Sh  =  Shape  factor,  3.32  for  sphere,  2.0  for  infinite  cylinder,  0.88  for  slab,  2.5  for 
cylinder  with  length  approximately  equal  to  diameter 

The  main  unknowns  in  the  F-K  analysis  are  the  thermal  conductivity,  the  heat  of  reaction,  and  the 
AiThenius  rate  constant  parameters,  A  and  Ea,  which  give  the  rate  constant  (k)  value  at  any 
temperature  as  k  =  A  exp(-Ea/RT).  The  density,  radius,  and  shape  factor  are  usually  known.  If  all 
these  parameters  are  known,  accurate  predictions  of  critical  temperatures  can  often  be  made  (when  the 
situation  fits  the  assumptions  of  the  F-K  model). 

Conversely,  if  the  critical  temperature  is  experimentally  measured  and  the  thermal  conductivity  is 
known,  combinations  of  A  and  Ea  can  be  determined  which  satisfy  the  F-K  formula.  It  is  found  that 
if  Ea  is  estimated  and  the  A-factor  is  determined,  the  resulting  value  of  k  is  only  slightly  dependent 
upon  the  value  of  Ea  used.  Thus  we  can  determine  a  value  for  k  at  the  temperature  Tc  by  measuring 
the  critical  temperature  of  an  explosive  in  a  given  sample  size  and  shape. 

Odier  rate  constant  values  at  different  temperatures  can  be  determined  by  measuring  critical 
temperatures  for  other  sample  sizes.  A  fit  to  these  k(Tc)  data  points  vs.  1/T  yields  a  global  rate 
coefficient  which  will  describe  cook-off  behavior  for  this  explosive  over  a  wide  range  of  conditions. 

2.3  SEMENOV  MODEL. 

For  containers  of  hot  stirred  liquids,  the  temperature  profile  from  the  edge  to  the  center  of  the  material 
will  not  show  the  expected  variation  appropriate  to  the  F-K  model  of  heat  losses  only  by  conduction. 
The  situation  of  a  well  stirred  material  with  an  almost  constant  temperature  profile  across  the  material 
is  described  by  the  Semenov  model6  according  to  the  equation 

V  Q  p  A  exp(-Ea/(R  Tc))  =  a  SA  Tc2  /Ea  , 

where  V  =  Volume  of  explosive,  cm3 

Q  =  Heat  of  decomposition  reaction,  cal/g 

p  =  Density,  g/cm3 

A  =  Preexponential  factor,  s'1 

Ea  =  Activation  energy,  cal/mol 

R  =  Gas  constant,  cal/K/mol 

Tc  =  Critical  temperature,  K 

a  =  Heat  transfer  coefficient  of  wall,  cal/(cm2  s  K) 

SA  =  Surface  Area,  cm2 

In  this  model,  the  same  values  for  the  kinetic  factors,  density,  and  heat  of  decomposition  would  be 
used  as  in  the  F-K  model.  The  surface  area  and  volume  would  be  characteristic  of  a  specific 
explosive  sample  and  would  be  known.  The  remaining  unknown  is  the  overall  heat  transfer 
coefficient,  a.  The  work  of  Rogers,  Guiana,  and  Loverro7  discuss  values  for  a  ranging  from  0.0085 
for  an  aluminum  wall  to  0.0022  for  a  steel  wall. 
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2.4  GENERAL  PROPERTIES  OF  AN  AND  ANFO8. 

Crystalline  ammonium  nitrate  absorbs  water  from  the  air  and  forms  hard  cakes  which  makes  it  difficult 
to  handle  (e.g.,  spreading  as  a  fertilizer).  These  handling  difficulties  led  to  the  manufacture  of  prilled 
AN  (i.e.,  small— approximately  2  mm  diameter— balls  of  AN  created  by  spraying  AN  solution  in  a 

drying  tower).  After  about  1945,  the  residual  moisture  content  of  the  prills  was  allowed  to  increase, 
producing  a  liigher  porosity  AN  denoted  as  fuel-grade  ammonium  nitrate  (FGAN).  The  previously 
used  low  porosity  AN  was  called  explosive-grade  ammonium  nitrate.  The  available  porosity  in  FGAN 
of  about  0.07  cm3/g  is  almost  exactly  the  amount  needed  to  absorb  the  quantity  of  oil  to  yield  an 
oxygen  balanced  (stoichiometrically  balanced  to  CO;  and  H-,0  reaction  products)  ANFO  product.  The 
AN  prills  also  contain  unavailable  porosity  (i.e.,  not  available  to  absorb  oil)  which  is  also  important  in 
giving  the  ANFO  a  low  final  density  and  in  sensitizing  it  to  shock  initiation. 

AN  has  several  solid  state  phases  over  the  temperature  range  of  interest.  This  complicates  its  use  in 
applications  where  expansion  may  cause  difficulties.  These  phases  and  some  of  their  characteristics 
are  summarized  in  Figure  2-1.  Of  special  interest  are  the  expansions  and  heats  of  transition  that  occur 
at  the  IV  to  III  transition  at  32.3°C,  the  III  to  II  transition  at  84.2°C,  and  the  II  to  I  transition  at 
125. 2°C  before  melting  at  about  169°C. 

The  explosive  characteristics  of  FGAN  became  widely  recognized  alter  the  explosion  of  a  shipload  of 
over  3,000  tons  at  Texas  City  in  1947.  The  specific  product  that  exploded  was  FGAN  coated  with 
about  1  percent  wax  for  anti-caking  purposes.  (Actually  the  "wax”  was  a  mixture  of  paraffin,  rosin, 
and  petrolatum.)  However,  this  mixture  is  more  sensitive  than  either  pure  AN  or  the  oxygen  balanced 
mixture  with  5.5  percent  wax.  The  explosion  was  caused  by  several  factors:  the  relatively  high 
sensitivity  of  the  wax-coated  FGAN;  the  large  charge  size;  the  high  initial  temperature,  about  32°C; 
and  the  paper  bags  in  which  the  material  was  shipped.  It  is  likely  that  the  paper-AN  combination 
exhibited  a  lower  thermal  stability  than  the  FGAN  itself  and  it  self-heated,  possibly  starting  the  fires 
that  ultimately  lead  to  detonation  of  the  whole  shipload  of  FGAN.  Following  the  Texas  City  incident, 
the  anti-caking  wax  coating  was  replaced  by  a  coating  of  kieseiguhr. 

2.5  DECOMPOSITION  KINETICS  OF  AN/FUEL  MIXTURES. 

Although  the  literature  on  the  thermal  stability  of  AN  is  extensive,  few  studies  have  been  performed 
on  ANFO.  A  review  of  the  studies  of  AN/hydrocarbon  mixtures  is  given  below  with  a  brief  overview 
of  the  work  on  AN. 

Hainer9  studied  the  thermal  decomposition  of  AN  and  3  percent  wax-coated  fertilizer  grade  AN 
(FGAN)  using  weight  loss  and  N20  evolution  measurements  over  the  temperature  range  220-290°C  for 
AN  and  70-230°C  for  FGAN.  At  the  lower  temperatures,  samples  as  large  as  100  pounds  were  used 
for  gas  evolution  measurements,  resulting  in  rate  constant  determinations  as  small  as  3  x  10'9s  ‘  at 
132°C  for  AN  and  2-10  x  10'1:  s'1  at  70°C  for  FGAN.  No  other  measurements  at  these  low  reaction 
rates  are  available  for  either  material. 

Hainer’s  kinetic  results  (Ea  in  keal/mol  units)  for  FGAN  are 

k  =  3.7  x  10’°  exp(-34.8/RT)  s  '. 


4 


Figure  2-1.  Ammonium  nitrate  phase  characteristics. 


Hainer  concluded  that  die  reaction  rates  of  AN  and  FGAN  were  similar  at  high  temperatures  (>200°C), 
but  that  of  FGAN  was  much  larger  at  low  temperatures  (by  about  100  times  at  100°C).  However,  we 
noted  that  the  rate  constant  for  AN  at  low  temperatures  that  Hainer  used  for  comparison  was 
extrapolated  over  many  decades  using  an  uncertain  activation  energy. 

Oxley,  Kaushik.  and  Gilson10  studied  AN,  emulsion,  and  AN  mixed  with  5  wt%  mineral  oil  (denoted 
ANMO  to  differentiate  it  from  commercial  ANFO)  by  isothermally  decomposing  the  test  material 
sealed  in  glass  tubes  and  by  using  DSC  techniques.  ANMO  reacted  slightly  faster  than  neat  AN  at 
temperatures  higher  than  about  270°C.  However,  at  lower  temperatures,  the  ANMO  reaction  was 
faster  than  AN  only  for  the  initial  part  (40  percent  at  270°C,  decreasing  to  14  percent  at  230°C)  of  the 
decomposition.  Kinetic  parameters  were  determined  for  both  the  more  rapid  initial  (denoted  "early") 
and  the  slower  final  (denoted  "late")  portions  of  the  total  reaction.  The  early  kinetics  of  ANMO 
decomposition  were  characterized  over  the  temperature  range  230-370°C  by 

k  =  2.4  x  101!  exp(-35.2/RT)  s'1  . 

2.6  OVERVIEW  OF  ANFO  DATA. 

The  kinetic  data  for  the  two  studies  of  AN/fuel  are  plotted  in  Arrhenius  form  in  Figure  2-2  for 
comparison.  The  activation  energies  of  the  two  expressions  are  similar,  with  the  two  fines  offset  from 
each  other  by  about  a  decade. 


2.7  DECOMPOSITION  OF  AN. 

Extensive  literature  exists  on  AN  decomposition,  including  information  on  kinetics,  mechanisms, 
catalyst,  and  product  composition.  King  and  Bauer11  at  Queens  University,  Ontario,  reviewed  the 
literature  on  AN  and  published  a  series  of  twelve  reports.  Number  four  of  this  series  covers  thermal 
decomposition  mechanisms  and  kinetics  up  to  1974.  Report  number  one  covers  accidents  with  AN 
and  number  nine  covers  the  deflagration  to  detonation  transition  in  molten  AN. 

Robertson12  studied  isothermal  AN  decomposition  by  pressure  measurements  over  the  temperature 
range  243-36 1°C.  This  work  extended  to  higher  temperatures  than  are  typical.  The  results  were  fit 
into  the  Arrhenius  expression 

k  =  6.3  x  1013  exp(-40.5/RT)  s'1, 

with  an  estimated  uncertainty  of  ±2.5  kcai/mol  in  the  activation  energy.  These  rate  constant  results  are 
higher  than  reported  in  most  other  work  done  in  this  temperature  range. 

Hainer’s  kinetic  results  for  AN  (220-290cC)  from  the  same  study9  as  discussed  above  were 

k  =  8.4  x  1018  exp(-55.6/RT)  s'*. 

It  should  be  noted  that  these  data  for  AN  did  not  cover  as  extensive  a  temperature  range  as  for  the 
FGAN  data. 

Cook  and  Abegg13  studied  AN  decomposition  using  real  time  weight  loss  measurements  over  the 
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Figure  2-2.  Arrhenius  graph  of  literature  data  on  AN/fuel  mixtures. 
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temperature  range  218-267°C.  First  order  behavior  was  observed  in  all  experiments,  with  no 
autocatalysis.  Their  rate  constant  was 

k  =  1.9  x  10,:  exp(-38.3/RT)  s’1  . 


Cook  and  Abegg's  expression  was  corrected  for  AN  vaporization  by  King  and  Bauer11  to  give 

k  =  4.9  x  1011  exp(-37.0/RT)  s'1  . 

Keenan  and  Dimitnades14  used  a  gas  evolution  method  to  study  the  decomposition  of  AN 
from  230-265°C.  Their  results  gave 

k  =  3.0  x  1016  exp(-49.5/RT)  s'1  . 

Rosser,  Inami,  and  Wise15  studied  AN  decomposition  from  225-275°C  using  N;0  evolution 
rates  in  a  flow  system.  By  seeding  the  inlet  gas  flow,  they  showed  clearly  that  NH3  and  H;0 
inhibit  AN  decomposition,  and  HN03  promotes  AN  decomposition;  thus  the  overall  rate  of 
decomposition  depends  on  the  gas  composition  in  contact  with  the  AN.  Their  nominal  rate 
constant  was 

k  =  1.8  x  1013  exp(-41.0/RT)  s1  . 

Brower,  Oxley,  and  Tewari16  calculated  rate  constants  for  AN  decomposition  in  sealed  glass 
tubes  from  total  gas  evolution  and  nitrogen  fraction  measurements  over  the  temperature  range 
200-380°C.  Their  data  were  fit  to  separate  high  and  low  temperature  expressions 

k  =  8.5  x  10M  exp(-46.3/RT)  s'1  (high  temperature),  and 
k  =  1.3  x  10  s  exp(-28.2/RT)  s'1  (low  temperature). 

2.8  OVERVIEW  OF  AN  DATA. 

The  kinetic  data  for  AN  are  plotted  in  Arrhenius  form  in  Figure  2-3  for  comparison. 

Although  all  of  these  rate  constant  expressions  are  derived  from  data  which  agree  fairly  well 
(within  a  factor  of  ten)  in  the  temperature  region  of  overlap,  they  were  almost  all  measured  at 
relatively  high  temperatures.  If  these  expressions  are  extrapolated  to  the  lower  temperatures, 
which  are  of  interest  for  storage  stability  considerations  for  large  quantities  of  AN  or  ANFO, 
the  disagreement  is  much  larger. 
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Figure  2-3.  Arrhenius  graph  of  literature  data  on  AN. 
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SECTION  3 


EXPERIMENTAL 


3.1  COOK-OFF  TESTS. 

The  critical  temperature,  Tc  (i.e.,  the  lowest  temperature  at  which  an  explosive  sample 
exhibits  thermal  runaway  from  chemical  self-heating)  was  determined  for  three  different 
sample  sizes  of  ANFO. 

3.2  HENKIN-MCGILL  TESTS. 

The  Henkin-McGill  test17  is  the  smallest  thermal  runaway  test.  It  has  the  advantage  of  being 
small  enough,  using  samples  of  less  than  0. 1  g,  to  be  run  in  the  laboratory  environment  using 
an  apparatus  similar  to  those  described  by  Rogers.18 

Ground  ANFO  samples  were  sealed  inside  standard  aluminum  blasting-cap  shells  using  a 
hand-operated  hydraulic  press  and  specially  made  hollow  aluminum  plugs.  A  two-step 
pressing  procedure  was  used  to  compress  the  sample  and  then  to  expand  the  wall  of  the  plug 
to  create  a  seal  with  the  shell.  This  gave  the  sample  a  reproducible  geometry  of  a  disc  (0.65 
cm  diameter  with  an  average  thickness  of  0.071  cm)  so  it  could  be  considered  to  be  an 
infinite  slab  for  analysis  (with  shape  factor  =  0.88).  A  micrometer  was  used  to  measure  the 
thickness  of  each  disc. 

The  temperature  of  the  Wood's  metal  bath  was  controlled  by  an  Omega  Engineering,  Inc. 
Model  920  proportional  controller.  The  actual  bath  temperature  was  measured  immediately 
prior  to  each  test  using  an  Omega  K-type  stainless  steel  sheathed  thermocouple  and  a  Fluke 
model  52  digital  thermometer  readout.  The  sample  was  lowered  into  the  hot  bath  in  a  remote 
operation  using  a  compressed  air-dnven  piston/cylinder.  The  time-to-explosion  was  measured 
using  a  digital  stopwatch  from  the  initial  time  when  the  sample  was  immersed  in  the  molten 
metal  bath  to  the  time,  if  any,  at  which  the  sample  holder  ruptured.  The  timing  uncertainty 
was  estimated  to  be  about  plus  or  minus  one  second.  To  find  the  critical  temperature,  the 
time-to-explosion  tests  were  repeated  at  various  temperatures  until  the  minimum  temperature 
was  found  at  which  thermal  runaway  and  shell  rupture  occurred. 

3.3  SSCB  TESTS. 

Pakulak  and  Cragin 19  developed  a  small  cook-off  test,  called  the  super  small  cook-off  bomb 
(SSCB)  test,  using  a  5/8  inch  (1.6  cm)  i.d.  tube  with  witness  plates  bolted  directly  to  the  open 
ends  of  the  tube  to  contain  the  explosive.  In  the  work  reported  here,  a  modified  SSCB 
apparatus,  based  loosely  on  the  Navy  test,  was  used.  Thick-walled  aluminum  tubing  with  2.5 
cm  o.d.  and  1.6  cm  i.d.  was  threaded  and  sealed  using  stainless  steel  pipe  caps  on  each  end. 

A  1/16  inch  (0.16  cm)  diameter  stainless  steel  sheathed  thermocouple  was  installed  through 
the  bottom  end  cap  using  a  Swagelok  fitting  to  measure  the  axial  sample  temperature.  A 
diagram  of  the  apparatus  is  shown  in  Figure  3-1.  The  assembled  tube  containing  the 
explosive  was  heated  by  contact  with  a  split  and  spring-loaded  4  inch  (10  cm)  diameter  brass 
platen  containing  four  embedded  250  Watt  electric  cartridge  heaters.  Heating  time  from 
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Figure  3-1.  Schematic  of  SSCB  cook-off  apparatus. 
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ambient  was  about  10  minutes.  The  tests  reported  here  were  run  vented  to  the  atmosphere. 

The  thermocouple  located  inside  the  test  vessel  and  the  two  thermocouples  located  m  the 
platen  were  connected  to  a  Hewlett-Packard  Vectra  personal  computer  using  an  Omega 
Engineering,  Inc.  White  Box  interface  card  and  software.  The  software  was  programmed  to 
record  the  temperatures  at  one-minute  intervals.  The  data  acquisition  board  was  checked  prior 
to  each  test  using  a  precision  voltage  source  thermocouple  calibrator.  The  channels  typically 
agreed  within  ±0.1  °C.  The  desired  platen  temperature  was  also  programmed  into  the  Omega 
software  which  then  controlled  a  solid-state  relay  connected  to  the  electric  heaters.  ASCII 
data  files  were  imported  into  a  Quattro  Pro  spreadsheet  for  data  analysis  and  plotting. 

3.4  1 -LITER  COOK-OFF  TESTS. 

Charges  of  about  0.8  kg  of  ANFO  were  tested  using  glass  round-bottomed  flasks  as  charge 
containers  using  three  0.16  cm  diameter  stainless  steel  sheathed  thermocouples  located  at  the 
sample  center,  mid-radius,  and  wall.  The  sample  was  heated  using  a  stirred  hot  oil  bath  using 
Dow  Coming  200  heat  transfer  fluid,  as  shown  in  Figure  3-2.  Another  thermocouple  was 
inserted  into  the  oil  bath  to  monitor  its  temperature.  A  two-bladed  stirrer  was  used  to  obtain 
uniform  oil  temperature.  The  same  data  acquisition  and  data  reduction  methods  used  for  the 
SSCB  tests  were  used  here. 

3.5  32-LITER  COOK-OFF  TESTS. 

Larger  quantities  of  ANFO  were  used  to  perform  thermal  runaway  tests  on  the  material  at 
temperatures  below  the  melting  point  (165-169°C).  To  minimize  fragments  should  the  charge 
detonate,  a  test  vessel  was  fabricated  based  on  a  32-liter,  heavy-walled  glass  battery  jar  (41 
cm  diameter,  3 1  cm  deep)  with  a  lightweight  aluminum  lid.  Silicone  rubber-insulated, 
flexible  electric  heating  bands  were  used  on  the  top,  bottom,  and  sides  of  the  vessel  for 
heating.  Insulating  blanket  material  was  also  wrapped  around,  over,  and  under  the  vessel. 

Three  0.16  cm  diameter  thermocouples  were  located  in  the  charge  as  shown  in  Figure  3-3. 

The  thermocouple  signals  were  connected  to  voltage-to-current  modules  (Omega  Engineering, 
OM3  modules)  using  eight-pair  K-type  thermocouple  wire  cable  to  transmit  the  signals  over 
approximately  210  meters  (600  feet)  of  shielded  multi-pair  cable  to  the  instrumentation 
bunker.  These  tests  utilized  the  same  type  of  computer  and  interface  board  that  were  used  for 
recording  the  data  and  controlling  the  heaters  for  the  SSCB  and  I -liter  tests.  The 
analog-to-digital  board  was  reprogrammed  from  direct  processing  of  thermocouple  voltage 
inputs  to  processing  of  20  miiliampere  current  loop  signals.  Remote  control  software 
(Close-Up)  and  high-speed  modems  (Hayes  Optima  9600)  were  installed  on  the  computer 
located  at  the  cook-off  site  (used  to  control  the  test  temperature  and  record  the  data),  as  well 
as  on  computers  used  for  remote  monitoring  of  the  test  progress. 

The  absolute  calibration  of  the  thermocouples  and  voltage-to-current  modules  were  checked 
using  a  voltage  source  calibrator  just  prior  to  beginning  the  test.  Slight  linear  offsets  were 
applied  to  the  various  thermocouple  signals  in  the  data  acquisition  program  to  bring  them  into 
agreement  within  ±0.5°C. 


12 


Figure  3-2.  Schematic  of  1 -liter  cook-off  apparatus. 
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A  plywood  box  was  placed  over  the  whole  apparatus  to  provide  protection  from  weather. 
These  tests  were  performed  at  the  remote  EMRTC  Fast/Slow  Cook-off  site  which  IS  rated  for 
large  explosive  charges  and  has  continuous  video  surveillance  for  security. 

3.6  TEST  MATERIAL. 

Commercial  ANFO  obtained  from  ICI  Explosives  in  Albuquerque  was  used  for  all  tests.  The 

ANFO  was  used  as  received  except  for  grinding  it  for  Henkin-McGiil  tests  (where  the  prills 
were  too  large). 
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SECTION  4 


RESULTS 

4.1  DENSITIES  AND  MELTING  POINT. 

The  density  of  the  ANFO  was  measured  four  times  using  a  cup  of  known  volume  and  a 
precision  balance.  The  results  were  an  ANFO  density  of  0.85±0.01  g/cmJ.  The  average 
density  of  the  Henkin-McGill  samples  which  were  pressed  into  the  shells  for  testing  was 
measured  in  place  as  about  1.65  g/cm3. 

The  melting  temperature  of  this  ANFO  was  measured  using  about  2  grams  of  material  inside 
a  test  tube  heated  by  an  oil  bath.  Some  melting  (wetting)  was  observed  at  160°C,  but  melting 
was  not  complete  until  about  170°C.  Melting  tests  were  also  performed  using  a  Mettler  hot 
stage  with  microscopic  observation.  No  melting  was  observed  below  167.6°C,  with  flooding 
observed  by  169.6°C.  The  melting  point  is  thus  between  167  and  169°C. 

4.2  HENKIN-MCGILL  COOK-OFF  DATA. 

About  twenty-five  40±1  mg  samples  were  tested  to  determine  the  small-scale  critical 
temperature.  These  samples  had  an  average  thickness  of  0.071  cm.  The  experimental  data 
are  listed  in  Table  1  (Attachment  1)  and  shown  in  Figure  4-1  (go's  and  no-go's  are  denoted  by 
open  and  filled  symbols,  respectively).  Little  overlap  was  found  between  the  go's  and  no-go’s 
near  the  critical  temperature.  The  definition  of  the  critical  temperature  is  the  temperature  for 
the  lowest  go,  ignoring  that  there  may  be  some  no-go's  at  higher  temperatures.  Thus  the 
critical  temperature  for  this  sample  size  of  ANFO  was  determined  to  be  297UC.  The  positive 
tests  were  characterized  by  relatively  weak  events. 

4.3  SSCB  COOK-OFF  DATA. 

Two  SSCB  tests  were  run  on  approximately  25  g  samples  of  ANFO  in  1.6  cm  diameter 
aluminum  cylinders  at  temperatures  above  the  melting  temperature.  These  tests  were  tried 
with  the  hope  that  the  small  cross-section  of  this  vessel  would  prevent  convection  in  the 
liquid  ANFO.  The  results  of  the  first  SSCB  test  at  230°C  are  shown  in  Figure  4-2.  The  test 
temperature  was  reached  in  about  14  minutes.  After  about  37  minutes,  a  rise  in  sample 
temperature  is  shown,  followed  by  a  temperature  decrease.  Although  this  change  in 
temperature  appears  relatively  minor,  visual  monitoring  of  the  test  apparatus  via  a  video 
camera  showed  that  a  short  length  of  steel  tubing  connected  to  the  vent  fitting  of  the  vessel 
had  been  ejected.  Smoke  and  particles  escaped  from  the  vent  following  this  time,  indicating 
that  a  cook-off  had  occurred.  The  decrease  in  temperature  was  thought  to  be  caused  by 
ejection  of  a  significant  amount  of  the  liquid  ANFO,  exposing  the  tip  of  the  thermocouple  and 
leading  to  lower  temperature  measurements. 

The  second  SSCB  test  was  run  at  a  test  temperature  of  221°C,  as  shown  in  Figure  4-3  on  the 
following  page,  for  about  180  minutes,  followed  by  an  additional  1600  minutes  at  225CC.  No 
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Figure  4-1.  Henkin-McGill  data  for  40  mg  cook-off  tests  on  commercial  ANFO. 

The  filled  symbols  represent  no-go  results  and  the  open  symbols 
represent  go's.  The  critical  temperature  was  found  to  be  297°C. 
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Figure  4-2.  SSCB  cook-off  test  #1  for  ANFO.  An  event  occurred  at  about  38 
minutes  at  230°C. 
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Figure  4-3.  SSCB  cook-off  test  #2  for  ANFO.  No  thermal  runaway  was  observed 
at  225°C. 
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chemical  self-heating  was  observed. 


The  conclusion  reached  was  that  the  critical  temperature  for  this  sample  size  of  ANFO  was 
between  221  and  230°C,  about  225°C. 

4.4  1 -LITER  COOK-OFF  DATA. 

One-liter  isothermal  cook-off  tests  were  run:o  at  160,  163,  180,  and  205°C  without  observing 
thermal  runaway.  However,  visual  monitoring  of  the  tests  showed  that  considerable  gas 
evolution  was  occurring  during  the  test,  with  resulting  vigorous  stirring  of  the  molten  ANFO 
at  temperatures  above  170°C.  This  prevented  thermal  runaway  at  the  expected  temperature 
(approximately  180-190°C). 

The  Frank-Kamenetskii  (F-K)  model  assumes  that  heat  loss  is  only  by  conduction  processes. 
The  enhanced  heat  transfer  from  stirring  the  explosive  changes  the  situation  and  makes  the 
F-K  model  no  longer  applicable.  The  experimental  cook-off  test  data  are  still  valid,  but 
predictions  and  analysis  using  the  F-K  model  are  not  applicable. 

Approximate  Semenov  calculations  were  performed  to  estimate  the  critical  temperature  for  a 
well-stirred  1 -liter  test  sample  of  ANFO.  These  calculations  require,  in  addition  to  parameters 
used  in  the  F-K  model,  a  value  for  the  overall  heat  transfer  coefficient  at  the  wall  of  the 
container.  This  value  is  not  precisely  known.  However,  for  the  expected  range  of  the  heat 
transfer  coefficient,  Semenov  calculations  yielded  critical  temperatures  for  molten  ANFO  in 
the  2I6-233°C  range.  Since  this  temperature  range  is  higher  than  where  any  of  the  cook-off 
tests  were  performed,  the  calculation  supports  the  hypothesis  that  stirring  the  sample  from  the 
bubbling  phase  prevented  runaway  under  our  experimental  conditions. 

Because  the  primary  interest  was  in  solid  ANFO  and  due  to  uncertainties  in  analyzing  the 
data,  additional  1 -liter  tests  were  not  performed  at  higher  temperatures. 

4.5  32-LITER  COOK-OFF  DATA. 

The  first  32-liter  cook-off  test  was  performed  using  about  25  kg  (55  pounds)  of  ANFO.  In 
comparison  with  the  smaller  tests,  the  distance  between  the  charge  and  the  data  acquisition 
was  much  greater  in  this  case,  about  210  meters  (600  feet).  This  required  electronic  line 
drivers  to  send  the  thermocouple  signals  over  long  cables.  The  cook-off  test  was  initiated  at  a 
wall  temperature  of  145°C.  At  about  4.1  days  into  this  first  large  ANFO  test,  the  apparatus 
suffered  a  near-hit  from  a  lightening  strike,  and  all  of  the  control  electronics  were  destroyed. 
When  power  was  returned  to  the  apparatus,  shorted  heater  relays  caused  a  fire.  The  ANFO 
burned  slowly  without  incident. 

The  second  32-liter  test  setup  was  identical  to  the  ill-fated  first  test.  The  initial  wall 
temperature  was  set  to  151°C.  The  ANFO,  being  a  good  thermal  insulator,  required  about 
four  days  to  reach  the  test  temperature  throughout  the  charge,  as  shown  on  the  temperature 
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profiles  in  Figure  4-4.  No  heat  generation  was  observed  during  this  period,  so  the  wall 
temperature  was  increased  to  158°C  at  a  test  time  of  5.2  days,  as  shown  in  more  detail  in 
Figure  4-5.  The  charge  was  held  under  these  conditions  for  about  2.5  days,  during  which  the 
ANFO  temperature  gradually  increased,  but  no  self-heating  was  observed.  The  charge 
temperature  again  increased,  this  time  to  162°C,  just  below  the  melting  temperature  of  ANFO. 
The  test  was  continued  for  another  2.1  days  (from  7.7  to  9.8  days)  without  incident.  The 
temperature  was  then  increased  to  165°C  for  about  0.8  day,  and  finally  set  to  175°C. 

The  final  two  days  of  the  test  are  shown  in  Figure  4-6  on  expanded  scales.  The  mid-radius 
thermocouple  signal  was  lost  at  about  9.2  days  into  the  test  (due  to  a  thermocouple 
voltage-to-current 

module  failure,  again  from  lightening  in  the  area).  Not  ail  of  the  details  of  the  temperature 
profiles  are  understood.  Upon  increasing  the  wail  temperature  from  165  to  175°C,  the  center 
temperature  signal  began  to  increase  and  then  showed  an  abrupt  rapid  increase  when  that 
temperature  was  168°C.  The  ANFO  temperature  near  the  wall  would  have  been  slightly 
hotter.  A  rapid  runaway  reaction  occurred  which  ignited  the  ANFO.  Examination  of  the  test 
vessel  showed  that  about  90  percent  of  the  ANFO  burned,  leaving  a  fused  layer  about  3-4  cm 
thick  in  the  bottom  of  the  vessel.  The  critical  temperature  for  this  test  was  concluded  to  be 
about  168±5°C. 
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ture  data  for  32-liter  cook-off  test  #2. 
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Figure  4-5.  32-liter  cook-off  data  on  expanded  scale.  The  mid-radius  thermocouple 
signal  was  lost  at  about  9.2  days. 
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Figure  4-6.  32-liter  cook-off  data  on  expanded  scale.  The  sample  exhibited  rapid 
self-heating  soon  after  the  temperature  was  increased  at  10.75  days. 
The  ANFO  burned. 
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SECTION  5 


DATA  ANALYSIS 

5.1  THERMAL  CONDUCTIVITY  OF  ANFO. 

An  accurate  value  for  the  thermal  conductivity  of  ANFO  was  needed  to  predict  and  analyze 
cook-off  data  and  to  make  safety  calculations  for  storing  large  quantities  of  ANFO  and  AN 
prills.  An  experimental  thermal  conductivity  value  was  obtained  by  analyzing  the  temperature 
profiles  during  the  first  150  minutes  of  one  of  the  1 -liter  cook-off  tests20.  The  ANFO  did  not 
begin  to  melt  until  after  this  time.  The  solid  phase  changes  at  32,  84,  and  125°C  were 
evident  in  the  temperature  profiles  (as  is  melting  at  later  times).  The  endothermicity  of  these 
solid  phase  changes  makes  the  thermal  conductivity  analysis  uncertain,  but  we  need  an 
"engineering"  value.  By  lumping  the  phase  changes  in  with  the  heating  processes,  we 
obtained  an  effective  thermal  conductivity  for  practical  use  over  this  temperature  range. 

The  1 -liter  heat-up  temperature  profiles  used  for  this  analysis  are  shown  in  Figure  5-1.  Also 
shown  are  temperature  profiles  calculated  using  Williams21  spreadsheet  method  and  proposed 
values  for  the  thermal  conductivity.  These  calculations  also  required  values  for  the  density 
and  specific  heat  of  ANFO.  The  density  was  measured  as  0.85  ±  0.01  g/cm3.  A  specific  heat 
value  from  Urbanski22  was  used  of  1.5  J/(g  K),  as  measured  for  AN(III)  at  about  35-45°C. 

The  value  for  the  thermal  conductivity  of  ANFO  was  varied  until  the  calculated  temperature 
profiles  for  the  1 -liter  test  compared  acceptably  with  the  experimental  data.  Because  of  the 
phase  changes,  the  fit  was  not  as  good  as  had  been  possible  in  previous  work  for  emulsion 
and  other  materials.  The  final  value  determined  for  the  thermal  conductivity  of  ANFO  was 
0.1 1  W/(m  K),  compared  to  that  of  high  density  AN  of  about  0.24  W/(m  K).23 
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Figure  5-1.  Thermal  conductivity  fit  to  1 -liter  cook-off  data  on  ANFO.  The 

effective  ANFO  thermal  conductivity  was  varied  to  fit  the  experimental 
data.  Calculations  were  performed  using  0.11  W/(m  K)  for  the  thermal 
conductivity,  1.5  J/g  for  the  specific  heat,  and  0.85  g/cm3  for  the 
density. 
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5.2 


KINETIC  ANALYSIS. 


The  three  experimental  critical  temperatures  from  the  Henkin,  SSCB,  and  32-liter  cook-off 
tests  were  analyzed  using  the  Frank-Kamenetskii  formula  to  give  reaction  rate  constants  at 
each  of  the  three  temperatures.  The  following  parameter  values  were  used  in  the  analysis: 

•  density  of  1.65  g/cm3  (measured)  for  Henkin  and  0.85  g/cm3  (measured)  for  SSCB  and 
32-liter 

•  thermal  conductivity  of  0.24  W/(m  K)  for  Henkin  and  0.11  W/(m  K)  for  SSCB  and  32-liter 

•  heat  of  reaction  of  1260  J/g  (300  cal/g) 

•  shape  factors  of  0.88  for  Henkin,  2.5  for  SSCB,  and  2.7  for  32-liter 


This  analysis  gave  the  following  results: 


Test 

Weight 

Radius 

Critical 

Temperatur 

e 

Shape 

Factor 

Rate 

Constant 

Henkin 

0.04  g 

0.035  cm 

297°C 

0.88 

1.2  x  10V 

SSCB 

25  g 

0.79  cm 

225°C 

2.5 

2.5  x  10V 

32-liter 

22.7  kg 

20.5  cm 

168°C 

2.7 

5.7  x  10V 

The  expression  fitting  these  points  (Ea  in  kcal/mol)  is  given  below  and  shown  in  Figure  5-2. 
ANFO  k  =  5.3  x  101S  exp  (-46. 04/RT)  s'1 


Also  shown  in  Figure  5-2  are  lines  representing  the  ANMO  data  of  Oxley,  Kaushik,  and 
Gilson  and  the  FGAN  data  of  Hainer  over  their  respective  temperature  ranges.  The  agreement 
is  excellent,  even  though  the  activation  energies  of  the  two  earlier  studies  are  lower  than  that 
of  the  current  work.  The  current  result  for  ANFO,  the  two  AN/fuei  results,  and  selected  rate 
constants  for  AN  are  compared  in  Figure  5-3.  The  agreement  between  our  Henkin  point  with 
the  line  from  Oxley  et  al.  is  excellent,  as  is  the  agreement  between  our  SSCB  point  and 
Hamer's  FGAN  data.  Our  ANFO  data  show  a  slower  rate  of  decomposition  than  did  Hamer's 
FGAN  at  about  165°C.  Assuming  that  ANFO  decomposes  faster  than  AN  by  a  factor  of  2-5 
times,  the  comparison  with  all  of  the  AN  work,  except  that  of  Robertson  and  Brower  et  al. 
(low  T),  is  also  excellent. 

5.3  TOPAZ  SIMULATIONS. 

To  test  our  ANFO  kinetic  and  thermal  parameters,  we  performed  simulations  of  the  SSCB  and 
32-liter  cook-off  tests  using  a  PC  version  of  the  thermal  code  TOPAZ.  TOPAZ-2D  is  a 
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Figure  5-2.  Arrhenius  graph  of  ANFO  decomposition  data.  The  filled  points  are  the 
experimental  data  from  this  work  The  expressions  of  Oxley  et  al  and 
of  Hamer  are  shown  for  comparison. 
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Figure  5-3.  Comparison  of  current  and  literature  data  on  ANFO  and  AN. 
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In  simulating  the  SSCB  test  at  230°C,  we  initially  used  the  ANFO  decomposition  rate 
constant  determined  in  this  work  (denoted  as  "standard").  The  TOPAZ  calculation  showed 
thermal  runaway  slightly  faster  than  observed  (35  minutes  vs.  40  minutes  observed). 

Although  this  is  excellent  agreement,  a  few  additional  runs  were  made  to  determine  the 
sensitivity  of  the  calculations  to  changes  in  the  rate  constant.  When  the  rate  constant  was 
divided  by  2,  no  runaway  was  computed.  When  it  was  reduced  by  25  percent,  runaway  was 
calculated  at  45  minutes.  Thus  the  rate  constant  expression  is  probably  accurate  within  about 
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10  percent  at  230°C. 


At  220°C,  where  no  SSCB  runaway  was  experimentally  observed,  the  standard  rate  constant 
expression  gave  no  computed  runaway. 

In  the  TOPAZ  simulations  of  the  32-liter  cook-off  tests,  the  standard  rate  constant  gave 
runaway  slightly  slower  than  was  observed.  Multiplying  and  dividing  the  rate  constant  by  2 
and  then  multiplying  by  5  gave  results  that  indicate  that  the  rate  constant  expression  may  be 
slightly  low,  by  about  a  factor  of  2,  in  the  temperature  region  of  150-165°C.  This  uncertainty 
is  about  the  same  magnitude  as  the  square  points  which  represent  the  data  in  Figure  5-2. 

Overall,  agreement  between  calculated  and  experimental  time  to  runaway  is  excellent  for  the 
SSCB  tests  and  good  for  the  32-liter  test.  The  temperature  profiles  from  the  32-liter  test 
indicate  that  the  experiment  is  more  complex  than  the  simple  model. 

5.4  APPLICATIONS  TO  LARGE  CHARGES. 

With  reason  to  think  that  extrapolation  of  the  ANFO  thermal  stability  model  calculations  on 
large  storage  quantities  should  be  reasonably  valid,  we  performed  some  calculations  on  field 
applications. 

5.5  F-K  CALCULATIONS. 

The  variation  of  critical  temperatures  for  ANFO  cylinders  and  spheres  with  increasing  radius 
was  calculated  using  the  simple  F-K  formula  and  the  input  parameters  discussed  above. 

Figure  5-4  shows  these  minimum  runaway  temperatures  for  cylinders.  The  critical 
temperatures  lie  in  the  range  from  160  to  140°C  for  cylinders  with  radius  from  about  25  to 
150  cm.  As  a  point  of  reference,  a  55-gallon  drum  has  a  radius  of  about  28  cm  and  would 
weigh  about  175  kg  (385  pounds)  when  filled  with  ANFO.  The  critical  temperature  for  such 
a  drum  would  be  158°C.  A  1-ton  cylinder  with  diameter  equal  to  length  would  have  a 
diameter  of  55  cm  and  a  critical  temperature  of  147°C. 

Critical  temperatures  similarly  calculated  using  the  F-K  formula  for  spherical  quantities  of 
ANFO  are  shown  vs.  charge  weight  on  a  logarithmic  scale  in  Figure  5-5.  The  linear 
relationship  shown  in  this  figure  indicates  that  the  critical  temperatures  decrease  by  about 
10°C  for  every  ten-fold  increase  in  the  mass  of  ANFO.  Note  that  these  calculations  were 
performed  for  spherical  configurations  of  ANFO,  but  that  the  results  are  not  very  sensitive  to 
the  shape  of  the  mass,  as  long  as  it  is  relatively  compact.  The  critical  temperature  for  2,000 
tons  is  seen  in  Figure  5-5  to  be  about  110UC  (230°F). 

Specific  calculations  were  performed  for  hemispheres  of  2,400  and  4,800  tons  of  ANFO. 

These  charges  would  have  radii  of  1070  and  1347  cm,  respectively,  at  a  density  of  0.85 
g/cm3.  F-K  calculations  were  performed  using  these  radii,  the  model  parameters  given  above, 
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and  a  shape  factor  of  3.33  (sphere),  ignoring  that  the  charges  were  only  half  spheres.  This 
assumes  that  the  interface  between  the  ANFO  and  the  ground  is  approximately  adiabatic  (no 
heat  loss),  which  is  probably  a  good  assumption.  Critical  temperatures  are  given  below  for 
the  standard  rate  constant,  and  for  the  rate  constant  multiplied  and  divided  by  a  factor  of  two. 
Since  the  uncertainty  in  the  rate  constant  is  thought  to  be  less  than  this  factor,  a  conservative 
range  of  critical  temperatures  were  calculated  to  include  this  and  other  uncertainties. 


Critical  Temperature 


ANFO  Charge  Size 

Standard 

Std.  x  2 

Std.  /  2 

2,400  tons 

107°C 

103°C 

1 12°C 

4,800  tons 

104°C 

100°C 

109°C 

Thus  the  final  F-K  results  are  107±4°C  and  104±4°C  for  the  2,400  and  4,800-ton  charge  sizes. 

5.6  TOPAZ  CALCULATIONS. 

A  TOPAZ  model  of  a  2,400-ton  hemisphere  of  ANFO  (radius  =  1067  cm)  was  created  using 
42  nodal  points  defining  31  volume  elements.  Because  of  the  large  size  of  the  charge,  the 
thermal  effects  of  the  container  were  ignored.  (Fiberglass  is  a  good  heat  conductor  compared 
with  ANFO.) 

The  first  calculations  were  made  using  a  constant  elevated  temperature  boundary  condition  at 
the  outside  surface  of  the  charge.  Runs  were  made  with  this  temperature  set  to  110,  125,  and 
MOT,  with  the  initial  ANFO  temperature  at  35°C.  The  calculations  were  run  to  simulate  a 
30-day  period  under  these  conditions.  No  thermal  runaway  occurred  in  these  calculations;  in 
fact,  the  ANFO  is  such  a  good  thermal  insulator  that  the  temperature  increase  was  limited  to 
the  material  near  the  surface  of  the  container.  Another  simuiadon  was  run  with  the  surface 
temperature  at  165°C,  well  above  the  critical  temperature,  and  the  inidal  ANFO  temperature  at 
80°C.  Although  the  temperature  of  the  material  near  the  surface  exceeded  100°C,  no 
significant  heat  generadon  or  runaway  was  seen. 

We  concluded  that  heating  in  this  manner  (i.e.,  by  conduction  from  the  outside  surface  of  the 
hemisphere)  cannot  create  a  runaway  situation  in  any  credible  scenario. 

5.7  ANFO  FIRE  HAZARDS. 

Another  possible  hazard  situation  would  be  a  fire  in  part  of  the  ANFO  charge  which  might 
heat  the  adjacent  ANFO  to  thermal  runaway  or  accelerate  to  a  deflagration  to  detonation 
(DDT)  transition.  Wang:o  previously  performed  some  pressurized  burning  tests  on  ANFO  and 
ANFO  with  10  wt%  aluminum  (Alcoa  1620A).  The  material  was  ignited  using  Reynolds 
Industries  Systems,  Inc.  SQ-80  thermite-filled  exploding  bridgewire  ignitors  located 


approximately  2.5  cm  above  the  samples.  These  tests  were  considered  preliminary,  but 
additional  tests  were  not  performed. 


"No-Go"  Pressure 

"Go"  Pressure 

(MPa) 

(MPa) 

ANFO 

4.1 

4.8 

ANFO  +  A1 

1.73 

2.76 

Note,  however,  that  in  a  later  test  on  ANFO  at  5.5  MPa,  the  combustion  failed  to  propagate 
through  the  complete  sample. 

A  few  linear  bum  rate  (LBR)  tests  were  performed  on  ANFO  +  10  wt%  aluminum  using 
squib  ignitors,  2.5  cm  diameter  glass  tubes,  and  fuse  wires  for  measuring  the  burning  front 
arrival  at  selected  locations  along  the  tube.  These  data  should  also  be  considered  preliminary. 
Tests  were  performed  at  two  pressures. 


Initial  Pressure 

Average  Pressure 

LBR 

(MPa) 

(MPa) 

(cm/s) 

5.52 

7.41 

1.626 

4.14 

5.17 

0.798 

These  data  show  two  main  points.  First,  for  small  samples,  ANFO  at  room  temperature  and 
atmospheric  pressure  does  not  readily  bum.  This  may  not  be  true  if  scaled  to  large  quantities 
of  material  where  heat  losses  would  be  reduced.  In  our  32-liter  cook-off  tests,  we  burned 
most  of  the  ANFO,  but  the  material  was  initially  quite  hot,  just  below  the  melting 
temperature.  Once  ignited,  melting  of  ANFO  at  approximately  169°C  also  may  help  since  the 
molten  material  may  flow  away  from  the  fire  area.  All  of  the  ANFO  was  not  consumed  in 
the  fires  that  resulted  from  burning  the  32-liter  test  charges. 

The  second  point  from  Wang's  data  is  that,  at  least  for  the  aluminized  ANFO  that  was  tested, 
the  bum  rate  is  rather  high  at  1-2  cm/s  under  the  high  pressure  conditions  tested.  However, 
in  ambient  fires,  these  pressures  would  be  difficult  to  achieve. 


35 


SECTION  6 


CONCLUSIONS 

Overall,  the  work  reported  here  showed  ANFO  to  be  a  highly  thermally  stable  energetic 
material,  suitable  for  use  in  large  conventional  explosive  charges.  The  data  measured  in  this 
work  allow  predictions  of  minimum  thermal  runaway  conditions  for  ANFO  charges  of  various 
sizes  and  for  various  initial  conditions.  If  the  initial  ANFO  temperature  loaded  into  a 
container  such  as  a  hemisphere  of  2,400  tons  (10  meters  radius,  70  feet  in  diameter)  is  near 
ambient,  our  model  indicates  that  no  credible  external  heating  is  likely  to  bring  the  material  to 
a  runaway  condition.  The  mam  remaining  hazards  are  from  a  large  fire  at  the  surface  or 
inside  the  ANFO,  or  in  some  other  nearby  fuel,  or  from  accidental  detonation  of  some  other 
explosive,  such  as  a  high  explosive  booster  charge,  near  the  ANFO. 

Specifically,  the  following  tasks  were  accomplished  in  this  program: 

•  The  effective  thermal  conductivity  of  solid  ANFO  was  measured  as  0.11  W/(m  K) 
using  a  specific  heat  value  of  1.5  J/(g  K). 

•  Isothermal  cook-off  tests  were  performed  on  samples  ranging  from  0.04  g  to 
approximately  25  kg. 

•  These  cook-off  data  have  been  used  to  give  the  thermal  stability  of  ANFO  from  168  to 
297°C — a  range  covering  over  five  decades  in  reaction  rates.  The  resulting  Arrhenius 
rate  constant  expression  was: 

k  =  5.3  x  1015  exp(-46.04/RT)  s'1  . 

These  data  compare  favorably  with  previous  studies  on  AN/fuel  mixtures. 

Minimum  thermal  runaway  temperatures  were  calculated  for  large  storage  quantities  of 
ANFO. 
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Appendix 


Henkin  Data  for  ANFO 
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ATTN:  TECH  LIBRARY 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  GBINNINGER 

SRI  INTERNATIONAL 

ATTN:  DR  JIM  GRAN 
ATTN:  EUTHE 
ATTN:  J  GIOVANOLA 
ATTN:  J  SIMONS 
ATTN:  M  SANAI 

TECH  REPS,  INC 

ATTN:  F  MCMULLAN 

TECHNICO  SOUTHWEST  INC 
ATTN:  S  LEVIN 

TITAN  CORPORATION 
ATTN:  J  ROCCO 
ATTN:  J  THOMSEN 


TITAN  CORPORATION  (THE) 
ATTN:  D  HATFIELD 
2  CY  ATTN:  F SAUER 
ATTN:  LIBRARY 

TRW  INC 

ATTN:  TIC 

TRW  S  I  G 

ATTN:  N  GUILES 

TRW  SPACE  &  DEFENSE  SECTOR 
ATTN:  W  WAMPLER 

W  J  SCHAFER  ASSOCIATES,  INC 
ATTN:  DYOUMANS 
ATTN:  W  BUITENHUYS 

WASHINGTON  STATE  UNIVERSITY 
ATTN:  PROF  Y  GUPTA 

WEIDLINGER  ASSOC,  INC 
ATTN:  H  LEVINE 

WEIDLINGER  ASSOCIATES.  INC 
ATTN:  TDEEVY 

WEIDLINGER  ASSOCIATES.  INC 
ATTN:  I  SANDLER 
ATTN:  M  BARON 

WILLIAMSON  AIRCRAFT  CO 
ATTN:  WILLIAMSON 
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